A cold-atom quantum simulator for SU(2) Yang-Mills lattice gauge theory 
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Non-abelian gauge theories play an important role in the standard model of particle physics, and 
unfold a partially unexplored world of exciting physical phenomena. In this letter, we suggest 
a realization of a non-abelian lattice gauge theory - SU(2) Yang-Mills in 1 + 1 dimensions, using 
ultracold atoms. Remarkably, and in contrast to previous proposals, in our model gauge invariance 
is a direct consequence of angular momentum conservation and thus is fundamental and robust. 
Our proposal may serve as well as a starting point for higher dimensional realizations. 



The importance of gauge theories, being in the heart of 
the standard model of modern particle physics, cannot be 
over estimated. Gauge fields give rise to the long range 
causal interactions between matter particles. However, 
while Quantum Electrodynamics (QED) is an abelian 
gauge theory, the strong interactions are described by 
Quantum Chromodynamics (QCD), a non-abelian SU (3) 
Yang-Mills gauge theory. This gives rise to many sig- 
nificant differences. For example, non-abelian theories 
involve the effect of quark confinement [T] [5], which is 
responsible for the forces that "bind quarks together", 
forbidding the existence of free quarks. This gives rise 
to the familiar structure of hadrons. This effect, as well 
as other non-perturbative phenomena of 3+1-d dimen- 
sional non-abelian models, required the development of 
new methods and techniques, such as lattice gauge theory 
[H [3HS] . It has been helpful to study fundamental QCD 
effects using simpler models that manifest the same es- 
sential ingredients. For example, confinement of quarks 
could be examined within the confining phase of (abelian) 
compact-QED [TJ [3J H3 [7] , which already in 2+1 dimen- 
sions gives rise to flux loops and plaquette interactions. 
Another approach is to use the 1+1-d version of SU(N C ) 
Yang-Mills theories (such as QCD 2 , the 1+1-d version 
of QCD); 1+1 dimensional models have been used in nu- 
merous nonperturbative methods to study the hadronic 
spectrum of such gauge theories [SHU] • 

The field of Quantum simulations [T5HT7] has been 
extensively developed theoretically and experimentally 
with the aim of advancing new quantum computational 
approaches to many body systems, in particular in the 
context of condensed matter physics. Recently, it has 
been realized that quantum simulations could be used 
also for exploring High Energy Physics (HEP) models 
and effects 18 . Simulations of HEP models are in gen- 
eral more demanding, as compared with condensed mat- 
ter ones for several reasons. Matter and gauge fields are 
described by fermionic and bosonic fields, and thus the 
use of several atomic species is needed. Furthermore, 
the continuum limit of the models must include rela- 
tivistic symmetries, which can be simulated with non- 



relativistic atoms by using lattices, as in lattice gauge 
theories. There have been several suggestions for the 
simulation of rclativisitic field theories involving bosons 
[PW421) . fermions, free or interacting with external (clas- 
sical) gauge fields [2"2"H2"4"] . or fermions interacting with 
(bosonic quantum) gauge fields in 1+1-d [25] . 

The simulation of dynamic gauge theories is even more 
challenging. First, gauge invariance must be respected. 
Moreover, it requires a special form of the many-body 
interactions, which are usually not available in most sys- 
tems. In particular, lattice models involve plaquette in- 
teractions among the links of the lattice. Second, Gauss 
law (involving gauge bosons and fermions) has to be im- 
plemented as a constraint. The special interactions could 
be obtained as a low-energy effective gauge invariant the- 
ory of the original system. 

Very recently, realizations of abelian dynamic gauge 
theories, employing trapped atoms in optical lattices, 
have been proposed, simulating 2+1-d Kogut-Susskind 
(KS) abelian compact QED, that manifests confinement 
of charges, using either BECs [21] or single atoms (a trun- 
cated version) [27] . These can also be extended to include 
dynamic matter, as proposed for the 1+1 dimensional 
Schwinger model [25] . which can be compared to exact 
available solutions, and for a 2+1 dimensional truncated 
KS model |25] . Simulations of other abelian gauge theo- 
ries [30ti32] have been proposed as well. 

In this letter we present a simulation scheme for an 
SU(2) Yang-Mills theory, where both the fermions (mat- 
ter) and gauge bosons are dynamical. The main idea is 
to introduce additional fermionic and bosonic fields (an- 
cillas) to obtain the desired non-abelian aspect. Remark- 
ably, and in contrast to previous proposals, in our model 
gauge invariance is fundamental and a direct consequence 
of angular momentum conservation, making this impor- 
tant property fundamental and robust. 

Lattice Yang-Mills Theory. The system we wish to 
simulate includes a non-abelian gauge field and dynamic 
fermionic matter. In ordinary lattice theory the gauge 
field degrees of freedom are defined on the lattice's links, 
whereas the matter fields are located on the vertices [31- 
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[S]. The gauge field is represented by unitary matrices U*, 
whose elements (U^) kl are constructed out of operators 
in the local (gauge field) Hilbert space. The index n 
labels the link (according to the vertex from which it 
emanates), and r the representation. In our SU(2) case, 
the fundamental representation is r = h and we shall 

1/2 

suppress the index r in this case. Thus U n = U n are 
2x2 matrices of operators. 

Non-ablian fields generally carry color charges, and 
hence, unlike in the abelian case, each link carries two 
different electric color fields, the left and right field char- 
acterized by the operators {L n a } , {R n , a }- Their differ- 
ence along a link can be interpreted as the color charge 
carried by it. These are, in fact, the left and right gen- 
erators of the group, and hence they must satisfy the 
(matrix) algebra of SU(2) [S3], 

[L a ,U r ]=\ r a U r ; [R a ,U r ] = U r \ r a (1) 

where {A^}, the representation matrices, are Aa -1 ^ 2 = 
\o a . As generators of SU{2), the left and right electric 
fields satisfy the algebra 

The left and right generators can be shown to commute 
with each other, and thus give rise to the same Casimir 
operator, L 2 = ^L a L a — R 2 . Thus, in SU(2) the gauge 

a 

field Hilbert space on a single link is characterized by 
three different integer quantum numbers, j,m,m', satis- 
fying 

L 2 \ jmm'} = R 2 \jmm!) = j (j + 1) \ jmm!) 
L z \ jmm'} = m \ jmm'} ; R z \ jmm'} = m! \ jmm!) (3) 

Finally, the r = 1/2 matter fields are introduced as 2- 
component spinors ip n defined at the vertices. The local 
color charges are defined by Q n . a = ^J^Yn k ( a a)ki Vv- 

k,l 

The local gauge invariance is manifested by the conserva- 
tion of Gauss's law at each vertex, L n a — R n -i. 
for each group index a separately. In local gauge trans- 
formations, one picks a group element V n for each ver- 
tex, and acts with it on the gauge and matter fields: 
ip n — > V n ip n ; U n — > V^U n V n+ i. This is a transformation 
in "group space", i.e. on group indices, and thus, in a 
gauge invariant Hamiltonian, all the group space indices 
(of matter and gauge fields) must be fully contracted (ef- 
fectively 'traced out'). Thus the simplest "pure gauge" 

terms are of the form cx Tr group {UiU 2 Ul Jjf\ , where the 
product is of group elements around a plaquette. Such 
terms give rise to the propagating effects in d + 1 di- 
mensions where d > 1, but are absent in 1+1-d, where 
the only possible interactions are gauge-matter ones. In 
the following we shall use the staggered fermions method 



in 1+1 dimensions 
Hamiltonian is 



where the gauge invariant 



H = J2 ( y L n + m(-l) n ^ n + i(3 &tU n i> n+1 - h.c.) 

n ^ 

(4) 

where g is the theory's coupling constant and m is the 
fermions mass. We shall denote by \vac) the zeroth order 
ground state in the strong coupling limit (g 2 3> /?, where 
the f3 part is treated as a perturbation). This state satis- 
fies \vac) = <g) 1 000) \iplipn = 1 - (-1)"). The 

links vertices n 

double filling in the odd vertices ("negative mass" ver- 
tices) is the " Dirac Sea" ; In the continuum limit of stag- 
gered fermions, two neighboring two-component spinors 
become a single four-component one. 

In order to simulate this Hamiltonian, one has to find 
an appropriate realization for the group elements {U n } 
and generators {£ n ,a} , {Rn.a} fulfilling the unitarity and 
algebra demands ( 1|2|3| ). This can be done using the 
Jordan-Schwinger map )37l I38j , connecting harmonic os- 
cillators (bosons) and angular momentum. Mapping be- 
tween SU(N) and bosonic systems [SHI HO] allows to ex- 
press the gauge field operators using bosonic atoms in the 
prepotential method [UJH2]. In this method, for SU(2), 
one defines four bosonic species on each link n: the first 
two are identified by the operators 01, ci2 on the left side, 
and the other two 61, 62 on the right side, constrained by 
Nl = Nr, where Nl = a\a\+a\a2 and Nr = b\bi+blb 2 . 
Then, one can express the unitary operators on each link 
as U = UlUr, where 



U L = 



VN L + 1 V 4 °i 



a\ -a 2 



U R = 



b\ bj\ 1 
-6a 61 J VNr + 1 
(5) 



Then, by identifying j = ^ = the generators 



k,l 



and L 2 = \ 



1 , R 2 



k.l 



Nr (Nr 
2 V 2 



l), equations 



fllpp! follow. 

Pure Gauge Simulation. First, we shall discuss the 
simulation of the gauge field, disregarding the fermions. 
The simulating system required for that is a set of opti- 
cal lattices |15j . Each minimum can contain two different 
(bosonic) atomic species out of four, Ai y2 or -Bi,2- The 
A, B minima alternate (see Fig. [lji,b). We assume that 
the energy levels of the bosonic modes on each minima 
are fairly separated, such that we can consider only the 
lowest one. Since nearest neighbor minima cannot con- 
tain similar atomic species, tunneling is eliminated, and 
thus the only remaining interactions are within the min- 
ima - scattering and number terms. Tuning the optical 
parameters and the chemical potential properly, one gets 
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vertex 
LinkH n+i 



X B,D 



FIG. 1. A schematic representation of the required bosons 
(a,b) and fermions (c,d) superlattice structure; (e) shows 
schematically the combination of all the species to form the 
links-vertices structure. We have exaggerated the distance 
between the wells in order to make the structure more clear. 



the Hamiltonian 



E 



Nr. 



Nr. 



N, 



R,n 



N, 



R,n 



(7) 

but since [He, Nl, u — ^V"i?.,n] = 0, if the constraint 
Nl,tl — N^ n is initially fulfilled at all links, we get the 
desired Hamiltonian, He = QfY^n- This corresponds 

n 

to a 1+1 dimensional SU(2) pure-gauge theory. How- 
ever, since this system has no dynamics at all, we would 
like to introduce some dynamic color charges (fermions). 

Dynamic Fermions Simulation. First, we show how to 
realize the fermionic mass term. In order to do that, we 
introduce two-component spinors ip n at every vertex (see 
Fig. [lj:), i.e. to the left of the A bosonic minima, with 
an alternating chemical potential, yielding the requested 
Hamiltonian H m = mj^(— 1)"^ ip n . This requires the 

n 

use of a superlattice, as indicated in Fig. [TJ;. Next, we 
introduce the nontrivial interaction term. In order to do 
that, we introduce ancillary fermionic species, "sitting" 
in the middle of the links, to the right of the A bosonic 
minima (see Fig. [TJi) . On each such "virtual vertex" we 
define a two-component spinor, Xm with the local Hamil- 
tonian H x = A^XnXn- Besides these spinors, we intro- 

n 

duce more bosonic species, and D\ t %, whose minima 
overlap with the ones of A 1)2 and B i 2 respectively (see 
Fig. [lji-b), serving as reference baths, all prepared in an 
identical coherent (i.e. BEC) state \a), where a G M, 
a > 1. 

The required boson-fermion interaction is obtained as 
a boson-assisted tunneling, where fermionic tunneling is 



accompanied by an internal boson change, specifically 

ff /=2i^EfW)i(^L(Xn)i + 

+ (xt) i (w R , n ) i .(<<P n+1 ) j +h.c.) (8) 



where 



W L = 



a\ci — a2c\ 
a\c2 a\c\ 



W R = 



_ ( bid, 



b\d 2 



—b 2 d\ b\d\ 



(9) 



The overlapping of the Wannier functions of the different 
atomic species (see Fig. [I]) naturally gives rise to such 
interactions. The minima of the same species are always 
separated by minima of other ones, and hence tunneling 
between neighboring same-species minima are unlikely. 
Thus, only the scattering processes contribute. 

In order to eliminate undesired processes and keep only 
the needed ones, one has to carefully select the hyper- 
fine levels of each of the atomic species such that only 
the desired processes would conserve angular momentum. 
Thus, gauge invariance will be a fundamental property of 
the system (and not an effective one), arising from the 
natural atomic angular momentum conservation. There 
are many possible choices for the magnetic levels to be 
used. If one, for example, chooses 7 Li atoms for the 
bosons and 40 K atoms for the fermions, there are many 
choices of uif values that should work (the selection rules 
we consider depend only on the magnetic level mp). 
A possible choice is m(A\) = 3,771(^2) = 2,m{B\) = 
-3,m(B 2 ) = -2,ro(Ci) = l,m(C 2 ) = -3,m(D x ) = 
-l,m(D 2 ) = 3,m(V>i) = 3/2, m^a) = -3/2,m(xi) = 
7/2, m(x2) = —7/2 [33]. A detailed discussion can be 
found in the supplementary part [44] . 

Another term which conserves angular momen- 
tum and arises in the 

lE^n (N^n-i + N L>n ), 



proposed setup is Hf = 
but this can be eliminated 



by a proper choice of the Hamiltonian parameters, 
as is explained later. Scattering terms of the form 
XnXn {Nl.u + ^Yr.ti) are also angular momentum con- 
serving and thus possible. However, they vanish in the 
relevant subspace, A4, where the dynamics takes place 
(see below). Finally, the only angular momentum allowed 
if) — x and x — X scattering processes vanish in A4 . The 
local — ip scattering leads to a constant (total number 
of fermions) in M.. 

Since Cj |a) = a|a), and c\\a) a\a) (because a E 
K, a 3> 1, and the same applies for dj), we get that 
effectively, within our subspace of interest, we can replace 
Wl, W r with aW L ,aW R , where W L = y/N L + 1U L and 
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Wr = Ur^Nr + 1. Thus, one effectively gets 

n,i,j \ 

+ (xL)i (U R , n ) i:j (^n+l), y/N R , n + 1 + (10) 

Note that this Hamiltonian is gauge-invariant, if we take 
into account the x s as well when doing a gauge transfor- 
mation. 

Next, assume that H x is the largest energy scale, i.e. 
A 3> <? 2 ,m, e. In that case, we can treat He + H m + 
Hf + Hf as perturbations to the large constraint H x . If 
we prepare initially the system with no x atoms at all, 
we adiabatically eliminate them and obtain an effective 
Hamiltonian expansion in this subspace [45] . which we 
denote by M.. In first order, we get He + H m + Hf. 
In second order, there are several possible contributions. 
Each virtual action of Hf requires a second operation of 
it on the same link, in order to return to M.. The first 
type of processes yield the Hamiltonian Hf = — Hf 
(see the supplemental part [H] for details). Choosing 
= 7, we eliminate both Hf, H'f as anticipated above. 

The other type of processes are of utmost interest. 
There, a fermion hops to the middle of a link, and 
then makes its way to the other side. Mathemati- 
cally (details in the supplemental part [44 ), one gets 

Hp = {^nV N L.n + lU ny /N R , n + + h.C.) 

n 

2 

where /3 = — This yields the full effective Hamil- 
tonian, to second order, H e ff = He + H m + Hp. The 
gauge- fermions coupling is folded within Hp. Now, let us 
perform on all the "real" fermions the canonical trans- 
formation ip n — ¥ i n ip n . This transformation leaves H m 
invariant, but introduces the missing complex factors to 

A single action of Hp on \vac) raises the flux on the 
relevant link to j = \{Nl = Nr = 1). Since a link in its 
ground state contains no A, B bosons, only the creation 
operators part of U, U+, would contribute, and one would 
get, before its operation, y/N + 1 \ vac) — \vac), and 
after its action, y/N + 1U \vac) = y/N + 1U+ \vac) = 
\/2U+\vac) = \/2U\ vac). Thus we conclude that 
Hp \vac) = if$Yli vPnUn'Pn+i ~~ h.c.) \vac) as desired. 

n 

One can also show (as done in the supplemental part 
[44] ) that a double operation of Hp on \vac) is equiva- 
lent to the simulated system's case as well, since the flux 
never exceeds j = 1/2 (if initially j < 1/2 everywhere. 
Furthermore, all the gauge invariant states (with j = 1/2 
charges) can be obtained by operations of the Hp terms 
on \vac). Thus, we conclude that if the system did not 
contain initially any any links with j > ~, this of course 
also applies for higher orders, which can be obtained by 
iterating Hp several times. Therefore, if the system does 



not contain initially any links with j > |, one can ef- 
fectively drop the square roots of number operators and 
write 

Hp = {^ n U n il> n+1 - h.c.) (11) 

n 

as in Q, and get an accurate simulation for all orders in 
Hp and all the parameters' regimes. 

Initial state preparation and possible measurements. 
Initially, one can prepare the system in the ground state 
of the strong limit, \vac), in which the C, D atoms are in 
the | a) state previously defined, and there are no other 
bosons all around the lattice. The fermions should be 
prepared as explained before. The fermion interactions 
should be switched off, i.e. /? = or e = 0, for example, 
by deepening the optical lattice minima. One can then 
create, using single addressing lasers HZ], charges 
with the appropriate flux tubes connecting them. One 
could either create mesons, whose length should be odd 
(since we are using staggered fermions and " quarks" and 
"anti-quarks" are on alternating vertices) or baryons, sit- 
ting on the same vertex, as described in |3S]. As long as 
no states with j > \ on any link are produced, the dy- 
namics would be exactly as in the simulated model. Then 
the fermionic dynamics can be turned on. Afterwards, 
one can change the parameters m, g, and if it is done adi- 
abatically, one can go from the strong to the weak limit 
and see the consequences (no phase transition is expected 
in 1 spatial dimension, of course) . Measurements can be 
done by probing the number of "real" fermions and A, B 
bosons over the lattice. Another possibility is to real- 
ize the Wilson-Loop area law measurement, proposed in 
[2U1 |4"B] for a non-abelian system. 

Before concluding, we shall emphasize again that in 
this model gauge invariance is fundamental and exact, 
inherited from angular momentum conservation. This 
makes the model robust against errors and corrections: 
as the symmetry is already manifested in the basic Hamil- 
tonian ([8]), the system cannot leave the gauge invariant 
subspace. 

Generalization to more dimensions. As mentioned be- 
fore, lattice gauge theories in d + 1 dimensions (d > 1) 
contain "magnetic" plaquette terms, unlike the 1 + 1 di- 
mensional case presented here. Thus the generalization 
to higher dimensions requires more complicated tech- 
niques |49] . Another important consequence of these 
terms is that for d > 1 the y/N + 1 operators encountered 
here cannot be avoided, since these closed flux loops do 
not allow remaining in the j ' < | subspace, or, phrased 
differently, since the flux is no longer confined in a line, 
it can spread and raise j. In this work we have consid- 
ered the simpler 1+1 dimensional model, although some 
of the ideas introduced here may be used as a basis to 
build higher dimensional versions. 

Finally, let us emphasize that it is fair to recognize that 
the simulation proposed here requires setups that are 
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much more complex than the ones required for the simu- 
lation of condensed matter models, something which will 
make the realization of the present proposal extremely 
challenging. However, the purpose of our work is to show 
that this, in principle, is possible, and to trigger the de- 
velopment of the required experimental techniques, since 
simulation of gauge theories may have a strong impact 
well beyond atomic and condensed matter physics. 

After the completion of this work, two other propos- 
als for quantum simulations of non-abelian gauge theo- 
ries with cold atoms have been suggested, for a strong- 
coupling Rishon-Link model [5U] and an SU(2) gauge 
magnet [ST] . 
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SUPPLEMENTAL MATERIAL 

Gauge invariance and angular momentum 
conservation 

Our model requires the creation of interactions of the 
form 

where 

W L = ( a j Cl - a2 f );Wn=( b h h h ) (13) 

as explained, such interaction term is obtainable if one 
takes advantage of hyperfme angular momentum conser- 
vation. The key point is that on each bose-fermi inter- 
action, one boson-fermion couple is created instead of 
another pair, which is annihilated. Angular momentum 
conservation (which is a kind of rotating-wave approxi- 
mation) forces the total mp of the created pair to equal 
the rriF of the annihilated one. Thus, one should specif- 
ically choose the hyperfme levels (mps) of the various 
atomic species such that only the desired processes can 
occur. 



This is done the following way. First, let us look at 
the interactions on the left side of each link, which in- 
volve four bosonic species (^1,2 5 £1,2) and four fermionic 
species (ipi t z, Xi.2)- The mp values should be picked such 
that: (1) The total m value of a pair created in one in- 
teraction will exactly equal the total m of the pair an- 
nihilated in this interaction, and not equal the total m 
in any other boson-fermion pair in the system; (2) The 
total m value of any pair not participating in the desired 
interactions does not equal the total m of any other pair, 
either participating in the interactions or not. 

After fulfilling this requirement, one has to examine the 
right side interactions. The fermionic quantum numbers 
are now inherited from the left side, but the four other 
bosonic species (.81,21 ^1,2) have to be seeked in the same 
manner, following the same two requirements. 

The entire process can be described using two summa- 
tion tables, for the left and right sides, with the fermions 
and bosons m values as the column and row headings. 
The required Hamiltonian terms determine which four 
pairs of elements of each of the tables must be equal, and 
all the other elements must be unique. Thus, due to an- 
gular momentum conservation, one gets only the required 
terms, and the non-species-changing scattering described 
by Hf in the paper, which is later cancelled by H',. As 
explained in the paper, most of the other terms do not 
contribute inside M. 

The following tables present the summation tables for 
the choice of ms given in the paper. The colors match the 
proper elements of the Wl, Wr matrices, given below. 



Left processes example 



w L = ( i ci - a *f ) 

V a 2 c 2 aic[ J 





m(Ai) = 3 


m(A 2 ) = 2 


mid) = 1 


m(C 2 ) = -3 


m(V>i) = 3/2 


9/2 


7/2 


5/2 


-3/2 


m(V> 2 ) = -3/2 


3/2 


1/2 


-1/2 


-9/2 


m( X i) = 7/2 


13/2 


11/2 


9/2 


1/2 


m( X 2) = -7/2 


-1/2 


-3/2 


-5/2 


-13/2 



Right processes example 



w R = ( b ' d K 6 ^ ) 

\ -b 2 4 hd\ ) 



(15) 
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m(Bi) = -3 


m(B 2 ) = -2 


m(Di) = -1 


m(D 2 ) = 3 


m(ipi) = 3/2 


-3/2 


-1/2 


1/2 


9/2 


m(ip 2 ) = -3/2 


-9/2 


-7/2 


-5/2 




m( X i) = 7/2 


1/2 


3/2 


5/2 


13/2 


m(xa) = "7/2 


-13/2 


-11/2 


-9/2 


-1/2 



The second order effective Hamiltonian terms 



Second order Hp actions 



As explained in the text, in the second order, there 
are several possible contributions to the effective Hamil- 
tonian. Each virtual operation of Hf requires a second 
operation of it on the same link, in order to return to M. . 

In the first type of processes, a "real" fermion ip hops 
to the left, over B bosons, to the middle of the link, and 
then jumps back: 



V2A 



E 



Here we shall elaborate on the second order action of 
Hp on | vac). The cases involving products of two differ- 
ent link operators are easily deduced to be equivalent to 
the real simulated theory, as in the single operation case. 
The ones which involve product of operators of the same 
link require a more careful consideration. First, consider 
the terms which involve U and on the same link: 



3 (^) fc V^L,n + 1 (U n ) kl y/N R , n + 1 (^n+l), \vac) (19) 



(Xt) t (U R ,n) lm y/N R , n + l(il; n+1 ) m j 

^(iV fl ,„-l + l)^ n (16) 



e 2 

V2> 



and similarly, a "real" fermion ip can hop to the right, 
over A bosons, and back, giving rise to a similar contri- 
bution. Together they add up to (disregarding a constant 
part, proportional to the total number of fermions in the 
system) 



H'f = -^$yi$n (^«,n-l + N L<n ) (17) 



Note that if by setting 7 



V2X 



in Hf , one gets Hf 



Hf = and both this unwanted Hamiltonians cancel 
each other and may be disregarded. 

The second type of processes is of more interest. There, 
a fermion hops to the middle of a link, and then makes 
its way to the other side. Mathematically, one gets 



V2X 



V2X 



n,i,j,l,m;\4>)^Ai \ 

I (xDt (u R , n ) lm Viv R , n + 1 (Vwi) m J + ^.c 



(18) 



but note that, as before, the rightmost viV~+T is just 1, 
and the AT+1 in the middle is 2. Moreover, (V>«.)j (Vvijfc = 

^'fc - (V4)fcfaMj- If n is Odd, (V>i) fc (V'n)j |«Oc) = 

(5jfc |uac), and we get zero, as in the simulated system. 
If n is even, ypX)k (V'n)i |«ac) = 0, and due to the uni- 
tarity of U n we get 



fi 2 \/N R , n + lip} l+1 ip n+1 \vac) = p 2 ipl +1 ip n+ i \vac) 



(20) 



as in the simulated system as well. 
Finally, consider terms of the form 



6 2 



setting j3 = — V one gets ff^ defined in the paper. 



2 W), + 1 V^fl,n + 1 (V-n+l), X 

W) fc VJVi.n + 1 (Un) kl y/Nk~n + 1 W>n+l), M (21) 

The three right square roots give us 2 as well, cancelling 
the denominator. The first (right) U operates on \vac), 
and thus it is U+. The second (left) one contributes 
U+ + U-. The U- contribution lowers the flux on the 
link to j = 0, and thus the leftmost y/N + 1 is 1 and we 
get an exact equivalence to the simulated system. The 
U+ contribution vanishes here as well as in the simu- 
lated system, due to canonical fermionic commutation 
relations. Thus, in conclusion, the second order of Hp is 
equivalent to the simulated system's case as well. Note, 
however, that all the gauge invariant states (with j = 1/2 
charges) can be obtained by operations of the Hp terms 
on | vac). Thus, we conclude that up to second order, Hp 
is equivalent to the simulated system, since it never pro- 
duces states with j > \, given that the system did not 
contain initially any any links with j > |. This of course 
also applies for higher orders, which can be obtained by 
iterating Hp several times. Therefore, if the system does 
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not contain initially any links with j > \, one can ef- write 
fectively drop the square roots of number operators and 

H = (V^n+i - h.c.) (22) 

n 

as in (4) in the letter, and get an accurate simulation for 
all orders in Hp and all the parameters' regimes. 



